Ultrasonic densitometer with opposed single transducer and transducer array

ABSTRACT

A scanning ultrasonic apparatus provides measurements of a bone at a number of spatially separated locations. These measurements may produce an image or may be used to automatically identify a measurement region of interest or to align a series of measurements with each other despite possible shifting in the bony member in between measurements.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a continuation-in-part of Ser. No. 09/094,073, filedJun. 9, 1998, now U.S. Pat. No. 6,027,449, which is acontinuation-in-part of Ser. No. 08/795,023, filed Feb. 4, 1997, nowU.S. Pat. No. 6,012,779, which is a continuation-in-part of Ser. No.08/466,495, filed Jun. 6, 1995, now U.S. Pat. No. 5,603,325, which is acontinuation-in-part of Ser. No. 08/397,027, filed Mar. 1, 1995, nowU.S. Pat. No. 5,483,965.

BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates to devices which are used for measuringthe density of members, such as bones, and more particularly to deviceswhich utilize ultrasonic acoustic signals to measure the physicalproperties and integrity of the members.

2. Description of the Prior Art

Various devices presently exist which may be used to measure thephysical properties and integrity of a member such as a bone.Non-invasive density measuring devices can be used to determinecumulative internal damage caused by micro-crushing and micro-fracturingoccurring in the bones of humans or animals such as race horses.Additionally, osteoporosis, or loss of bone mineralization, detection inhumans and its cure or prevention are increasingly becoming areas ofintense medical and biological interest. As the average age of the humanpopulation increases, a greater number of patients are developingcomplications due to rapid trabecular bone loss.

U.S. Pat. No. 3,847,141 to Hoop discloses a device for measuring thedensity of a bone structure, such as a finger bone or heel bone, tomonitor the calcium content thereof. The device includes a pair ofopposed spaced ultrasonic transducers which are held within a clampingdevice clamped on the bone being analyzed. A pulse generator is coupledto one of the transducers to generate an ultrasonic sound wave which isdirected through the bone to the other transducer. An electric circuitcouples the signals from the receive transducer back to the pulsegenerator for retriggering the pulse generator in response to thosesignals. The pulses therefore are produced at a frequency proportionalto the transit time that the ultrasonic wave takes to travel through thebone structure, which is directly proportional to the speed of the soundthrough the bone. The speed of sound through a bone has been found to beproportional to the density of the bone. Thus the frequency at which thepulse generator is retriggered is proportional to the density of thebone.

Another device and method for, establishing, in vivo the strength of abone is disclosed in U.S. Pat. Nos. 4,361,154 and 4,421,119 to Pratt,Jr. The device includes a launching transducer and a receivingtransducer which are connected by a graduated vernier and whichdetermine the speed of sound through the bone to determine its strength.The vernier is used to measure the total transit distance between thesurfaces of the two transducers.

Lees (Lees, S. (1986) Sonic Properties of Mineralized Tissue, TissueCharacterization With Ultrasound, CRC publication 2, pp. 207-226)discusses various studies involving attenuation and speed of soundmeasurements in both cortical and spongy (cancellous or trabecular)bone. The results of these studies reveal a linear relationship betweenthe wet sonic velocity and wet cortical density, and between the drysonic velocity and the dry cortical density. The transit times of anacoustic signal through a bone member therefore are proportional to thebone density. Langton. et al. (Langton, C. M., Palmer, S. D., andPorter, S. W., (1984), The Measurement of Broad Band UltrasonicAttenuation in Cancellous Bone, Eng. Med., 13, 89-91) published theresults of a study of ultrasonic attenuation versus frequency in the oscalcis (heel bone) that utilized through transmission techniques. Theseauthors suggested that attenuation differences observed in differentsubjects were due to changes in the mineral content of the os calcis.They also suggested that low frequency ultrasonic attenuation may be aparameter useful in the diagnosis of osteoporosis or as a predictor ofpossible fracture risk.

BRIEF SUMMARY OF THE INVENTION

The present invention provides an acoustic image of the human heel. Theimage can be used to provide greater insights into material andstructural variations within the heel, to locate a consistent region ofinterest on a given heel, or to develop a template which can be used toimprove the reproducibility of multiple measurements of a patient overseveral visits.

Specifically, the present invention provides an imaging ultrasonic bonedensitometer with at least one ultrasonic transducer arranged to measureacoustic signals modified by different portions of the bony member. Anelectronic data processor receives the electrical signals correspondingto the acoustic signals and processes the signals to determinecorresponding member variables related to the property of bony member atthe different locations. A display communicates with the data processorto provide a measure of the bony member at the positions. The membervariables may be attenuation, broad band ultrasonic attenuation (BUA),time of flight, speed of sound or a combination of these measurements.

It is thus one object of the invention to provide an ultrasonic bonedensitometer providing a spatially sensitive information about bonequality.

The display may be a graphic display providing an image of the bonymember, the image indicating the member variables as measured at thedifferent locations.

It is thus another object of the invention to provide a densiometricimage useful for evaluating bone quality.

The electronic data processor may operate to analyze the membervariables to identify a measurement region of interest in the bone. Themember variables within the region of interest may then be determined.

It is another object of the invention, therefore, to provide an imagingultrasonic bone densitometer where the image data can be used toaccurately locate a measurement region within the heel.

The densitometer may use an array of ultrasonic transducers providing afocused measurement of acoustic signals passing through a predeterminedlocation within the bony member. The electronic data processor may scanthe predetermined location through the bony member to provide a planaror volumetric image.

Thus it is another object of the invention to provide an ultrasonic bonedensitometer that may produce a high resolution densiometric image. Thepredetermined location may be shifted electronically to obtaininformation for a complete image both across the transmission path ofthe ultrasonic signals and at different depths within the bone along thetransmission path of the ultrasonic signal.

The electronic data processor may measure two locations within the bone,the first being within the trabecular region and the second at thecortical edge of the bony member.

Thus it is another object of the invention to provide a densitometerthat may make two spatially separate measurements indicating differenttypes of bone within its field of view.

The foregoing and other objects and advantages of the invention willappear from the following description. In this description, reference ismade to the accompanying drawings which form a part hereof and in whichthere is shown by way of illustration, a preferred embodiment of theinvention. Such embodiment does not necessarily represent the full scopeof the invention, however, and reference must be made therefore to theclaims for interpreting the scope of the invention.

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS

FIG. 1 is a perspective view of the ultrasound densitometer deviceconstructed in accordance with the present invention;

FIG. 2 is a perspective view of an acoustic coupler, two of which areshown in FIG. 1;

FIG. 3 is a front view of a transducer face from which acoustic signalsare transmitted or by which acoustic signals are received, the face ofthe other transducer being the mirror image thereof;

FIG. 4 is a schematic block diagram view of the circuitry of theultrasound densitometer device constructed in accordance with thepresent invention;

FIG. 5 illustrates the method of sampling a received waveform used bythe circuit of FIG. 4;

FIG. 6 is a schematic block diagram view of the circuitry of analternative embodiment of an ultrasound densitometer constructed inaccordance with the present invention;

FIG. 7 is a sample of an actual ultrasonic pulse and response from anultrasonic densitometer according to the present invention;

FIG. 8 is a sample plot of relative ultrasound pulse intensity overfrequency range;

FIG. 9 is a graph in frequency domain illustrating the shift inattenuation versus frequency characteristic of a measured object ascompared to a reference;

FIG. 10 is a perspective view of an alternative embodiment of thepresent invention showing a basin for receiving a patient's foot andhaving integral opposed ultrasonic transducers;

FIG. 11 is a plan view of a foot plate and toe peg used with theembodiment of FIG. 10;

FIG. 12 is a cross-sectional detail of the foot plate of FIG. 11 showingthe method of attaching the sliding toe peg of the foot plate;

FIG. 13 is a block diagram of a system for transporting the acousticcoupling liquid used in the embodiment of FIG. 10;

FIG. 14 is a schematic block diagram view of the circuitry of theembodiment of FIG. 10;

FIG. 15 is an exploded view of the underside of the foot basin of FIG.10 showing a c-clamp for holding the opposed ultrasonic transducers inprecise alignment and separation;

FIG. 16 is a perspective detailed view of the shank of the c-clamp ofFIG. 15 showing a lever for moving the separation of the transducersbetween an open and precisely separated closed position;

FIG. 17 is a cross-section of a human heel and ultrasonic transducers ofthe basin of FIG. 10 showing flexible liquid filled bladders surroundingthe transducers and providing a coupling path between the transducersand the heel;

FIG. 18 is a plot of the inverse of time of flight (TOF) for two boneconditions and broadband ultrasonic attenuation (BUA) as a function ofheel width showing their opposite functional dependencies;

FIG. 19 is a plot of bone quality versus bone width as might be obtainedfrom empirical measurement of multiple bone phantoms and as may be usedto eliminate bone width effects in the ultrasonic assessment of bonequality;

FIG. 20 is an exploded view of the elements of an ultrasonic detectorarray showing a driving mechanism for improving the resolution of theacquired data and the location of a piezoelectric film detector arrayabove a spatially offset connector;

FIG. 21 is a detailed perspective fragmentary view of the piezoelectricfilm detector with electrodes on its surface as communicating withconnector terminals via acoustically transparent conductors;

FIG. 22 is a detailed fragmentary view of the piezoelectric film showinga method of assembling the acoustically transparent conductors;

FIG. 23 is a detailed view of the face of the detector showing itsdisplacement by the driving mechanism of FIG. 20;

FIG. 24 is a figure similar to that of FIG. 17 showing use of thedetector array to provide focused reception at a point within apatient's heel;

FIG. 25 is a perspective view in phantom of a patient's heel showing araster scan pattern of a reception point within the heel to measurevolumetric bone density variations within an inner and outer portion ofthe os calcis;

FIG. 26 is a schematic representation of a data cube collected in thescanning shown in FIG. 25 with isodensity lines used to locate ameasurement region of interest;

FIG. 27 is a flow chart of the operation of the present invention inlocating a region of interest uniformly over several patient visits; and

FIG. 28 is a perspective view of an embodiment of the invention using afixed focus transducer array mechanically scanned to provide a pluralityof spatially separated measurements.

DETAILED DESCRIPTION OF THE INVENTION Caliper Embodiment

Referring more particularly to the drawings, wherein like numbers referto like parts, FIG. 1 shows a portable ultrasound densitometer 10 formeasuring the physical properties and integrity of a member, such as abone, in vivo. The densitometer 10 as shown in FIG. 1 includes a handle11 with actuator button 12. Extending linearly from the handle 11 is aconnection rod 13. The densitometer 10 also includes a fixed arm 15 andan adjustable arm 16. The fixed arm 15 preferably is formed continuouslywith the connection rod 13, and therefore is connected to an end 17 ofthe connection rod 13. The adjustable arm 16 is slidably mounted on theconnection rod 13 between the handle 11 and a digital display 18 mountedon the rod 13. The knob 19 may be turned so as to be locked or unlockedto allow the adjustable arm 16 to be slid along the connection rod 13 sothat the distance between the arms 15 and 16 may be adjusted.

Connected at the end of the fixed arm 15 is a first (left) transducer 21and at the end of the adjustable arm 16 is a second (right) transducer21. As shown in FIGS. 1 and 2, each of the transducers 21 has mounted onit a respective compliant acoustic coupler 23 to acoustically couple thetransducer to the object being tested. The acoustic coupler 23 includesa plastic ring 24 and attached pad 26 formed of urethane or othercompliant material. FIG. 3 shows a face 28 of the first (left)transducer 21 which is normally hidden behind the compliant pad 26 ofthe acoustic coupler 23. The transducer face 28 normally abuts againstthe inner surface 29 of the pad 26 shown in FIG. 2. The transducer face28 shown in FIG. 3 includes an array of twelve transducer elementslabeled a-l. The second (right) transducer 21 includes a face 28 whichis the mirror image of that shown in FIG. 3.

FIG. 4 generally shows in schematic fashion the electronic circuitry 31of the densitometer 10, which is physically contained in the housing ofthe digital display 18. An object 32 is placed between the twotransducers 21 so that acoustic signals may be transmitted through theobject. This object 32 represents a member, such as a bone, or somematerial with known acoustic properties such as distilled water or aneoprene reference block. As shown in the embodiment illustrated in FIG.4, the leftmost transducer 21 is a transmit transducer and the rightmosttransducer 21 a receive transducer. In fact though, either or both ofthe transducers 21 may be a transmit and/or receive transducer. Thetransmit and receive transducers 21 of the circuit of FIG. 4 areconnected by element select signals 36 and 37 to a microprocessor 38.The microprocessor 38 is programmed to determine which one of therespective pairs of transducer elements a through I are to betransmitting and receiving at any one time. This selection isaccomplished by the element select signal lines 36 and 37, which may beeither multiple signal lines or a serial data line to transmit theneeded selection data to the transducers 21. The microprocessor 38 isalso connected by a data and address bus 40 to the digital display 18, adigital signal processor 41, a sampling analog to digital converter 42,and a set of external timers 43. The microprocessor 38 has “on board”electrically programmable non-volatile random access memory (NVRAM) and,perhaps as well, conventional RAM memory, and controls the operations ofthe densitometer 10. The digital signal processor 41 has “on board”read-only memory (ROM) and performs many of the mathematical functionscarried out by the densitometer 10 under the control of themicroprocessor 38. The digital signal processor 41 specifically includesthe capability to perform discrete Fourier transforms, as iscommercially available in integrated circuit form presently, so as to beable to convert received waveform signals from the time domain to thefrequency domain. The microprocessor 38 and digital signal processor 41are interconnected also by the control signals 45 and 46 so that themicroprocessor 38 can maintain control over the operations of thedigital signal processor 41 and receive status information back.Together the microprocessor 38 and the digital signal processor 41control the electrical circuit 31 so that the densitometer 10 can carryout its operations, which will be discussed below. An auditory feedbackmechanism 48, such as an audio speaker, can be connected to themicroprocessor 38 through an output signal 49.

The external timer 43 provides a series of clock signals 51 and 52 tothe A/D converter 42 to provide time information to the A/D converter 42so that it will sample at timed intervals electrical signals which itreceives ultimately from the transmit transducer, in accordance with theprogram in the microprocessor 38 and the digital signal processor 41.The external timer 43 also creates a clock signal 53 connected to anexcitation amplifier 55 with digitally controllable gain. Timed pulsesare generated by the timer 43 and sent through the signal line 53 to theamplifier 55 to be amplified and directed to the transmit transducer 21through the signal line 56. The transmit transducer 21 converts theamplified pulse into an acoustic signal which is transmitted through theobject or material 32 to be received by the receive transducer 21 whichconverts the acoustic signal back to an electrical signal. Theelectrical signal is directed through output signal 57 to a receiveramplifier 59 which amplifies the electrical signal.

The excitation amplifier circuit 55 is preferably a digitallycontrollable circuit designed to create a pulsed output. Theamplification of the pulse can be digitally controlled in steps from oneto ninety-nine. In this way, the pulse can be repetitively increased inamplitude under digital control until a received pulse of appropriateamplitude is received at the receiver/amplifier circuit 59, where thegain is also digitally adjustable.

Connected to the receiver amplifier circuit 59 and integral therewith isa digitally controllable automatic gain control circuit which optimizesthe sensitivity of the receive transducer 21 and the amplifier circuit59 to received acoustic signals. The microprocessor 38 is connected tothe amplifier circuit and automatic gain control 59 through signal line60 to regulate the amplification of the amplifier circuit and gaincontrol 59. The amplified electric signals are directed through lead 61to the A/D converter 42 which samples those signals at timed intervals.The A/D converter 42 therefore in effect samples the received acousticsignals. As a series of substantially identical acoustic signals arereceived by the receive transducer 21, the A/D converter 42progressively samples an incremental portion of each successive signalwaveform. The microprocessor 38 is programmed so that those portions arecombined to form a digital composite waveform which is nearly identicalto a single waveform. This digitized waveform may be displayed on thedigital display 18, or processed for numerical analysis by the digitalsignal processor 41.

The densitometer constructed in accordance with FIGS. 1-4 can beoperated in one or more of several distinct methods to measure thephysical properties of the member, such as integrity or density. Thedifferent methods, as described in further detail below, depend both onthe software programming the operation of the microprocessor 34 as wellas the instructions given to the clinician as to how to use thedensitometer. The different methods of use may all be programmed into asingle unit, in which case a user-selectable switch may be provided toselect the mode of operation, or a given densitometer could beconstructed to be dedicated to a single mode of use. In any event, forthe method of use of the densitometer to measure the physical propertiesof a member to be fully understood, it is first necessary to understandthe internal operation of the densitometer itself.

In any of its methods of use, the densitometer is intended to be placedat some point in the process on the member whose properties are beingmeasured. This is done by placing the transducers 21 on the oppositesides of the member. To accomplish this, the knob 19 is loosened toallow the adjustable arm 16 to be moved so that the transducers 21 canbe placed on opposite sides of the member, such as the heel of a humanpatient. The outside surfaces of the pads 26 can be placed against theheel of the subject with an ultrasound gel 35 or other coupling materialplaced between the pads 26 and subject 32 to allow for improvedtransmission of the acoustic signals between the member 32 andtransducers 21. Once the transducers 21 are properly placed on themember, the knob 19 may be tightened to hold the adjustable arm 16 inplace, with the transducers 21 in spaced relation to each other with themember 32 therebetween. The actuator button 12 may then be pressed sothat acoustic signals will be transmitted through the member 32 to bereceived by the receive transducer 21. The electronic circuit of FIG. 4receives the electrical signals from the receive transducer 21, andsamples and processes these signals to obtain information on thephysical properties and integrity of the member 32 in vivo. Themicroprocessor 38 is programmed to indicate on the digital display 18when this information gathering process is complete. Alternatively, theinformation may be displayed on the digital display 18 when theinformation gathering process is completed. For example, the transittime of the acoustic signals through the member 32 could simply bedisplayed on the digital display 18.

Considering in detail the operation of the circuitry of FIG. 4, thegeneral concept is that the circuitry is designed to create anultrasonic pulse which travels from transmit transducer 21 through thesubject 32 and is then received by the receive transducer 21. Thecircuitry is designed to both determine the transit time of the pulsethrough the member 32, to ascertain the attenuation of the pulse throughthe member 32, and to be able to reconstruct a digital representation ofthe waveform of the pulse after it has passed through the member 32, sothat it may be analyzed to determine the attenuation at selectedfrequencies. To accomplish all of these objectives, the circuitry ofFIG. 4 operates under the control of the microprocessor 38. Themicroprocessor 38 selectively selects, through the element select signallines 36, a corresponding pair or a group of the elements a through l onthe face of each of the transducers 21. The corresponding elements oneach transducer are selected simultaneously while the remaining elementson the face of each transducer are inactive. With a given element, sayfor example element a selected, the microprocessor then causes theexternal timer 43 to emit a pulse on signal line 53 to the excitationamplifier circuit 55. The output of the excitation amplifier 55 travelsalong signal line 56 to element a of the transmit transducer 21, whichthereupon emits the ultrasonic pulse. The corresponding element a on thereceive transducer 21 receives the pulse and presents its output on thesignal line 57 to the amplifier circuit 59. What is desired as an outputof the A/D converter 42 is a digital representation of the analogwaveform which is the output of the single transducer element which hasbeen selected. Unfortunately, “real time” sampling A/D converters whichcan operate rapidly enough to sample a waveform at ultrasonicfrequencies are relatively expensive. Therefore it is preferred that theA/D converter 42 be an “equivalent time” sampling A/D converter. By“equivalent time” sampling, it is meant that the A/D converter 42samples the output of the transducer during a narrow time period afterany given ultrasonic pulse. The general concept is illustrated in FIG.5. The typical waveform of a single pulse received by the receivetransducer 21 and imposed on the signal line 57 is indicated by afunction “f”. The same pulse is repetitively received as an excitationpulse and is repetitively launched. The received pulse is sampled at asequence of time periods labeled t₀-t₁₀. In other words, rather thantrying to do a real-time analog to digital conversion of the signal f,the signal is sampled during individual fixed time periods t₀-t₁₀ afterthe transmit pulse is imposed, the analog value during each time periodis converted to a digital function, and that data is stored. Thus thetotal analog waveform response can be recreated from the individualdigital values created during each time period t, with the overallfidelity of the recreation of the waveform dependent on the number oftime periods t which are sampled. The sampling is not accomplishedduring a single real time pulse from the receive transducer 21. Instead,a series of pulses are emitted from the transmit transducer 21. Theexternal timer is constructed to provide signals to the sampling A/Dconverter 42 along signal lines 51 and 52 such that the analog valuesampled at time period to when the first pulse is applied to a giventransducer element, then at time t₁ during the second pulse, time t₂during the third pulse, etc. until all the time periods are sampled.Only after the complete waveform has been sampled for each element isthe next element, i.e. element b, selected. The output from the A/Dconverter 42 is provided both to the microprocessor 38 and to the signalprocessor 41. Thus the digital output values representing the complexwaveform f of FIG. 5 can be processed by the signal processor 41 afterthey are compiled for each transducer element. The waveform can then beanalyzed for time delay or attenuation for any given frequency componentwith respect to the characteristic of the transmitted ultrasonic pulse.The process is then repeated for the other elements until all elementshave been utilized to transmit a series of pulses sufficient to createdigital data representing the waveform which was received at the receivetransducer array 21. It is this data which may then be utilized in avariety of methods for determining the physical properties of themember. Depending on the manner in which the densitometer is beingutilized and the data being sought, the appropriate output can beprovided from either the microprocessor 38 or the signal processor 41through the digital display 18.

Because the ultrasonic pulsing and sampling can be performed so rapidly,at least in human terms, the process of creating a sampled ultrasonicreceived pulse can optionally be repeated several times to reduce noiseby signal averaging. If this option is to be implemented, the process ofrepetitively launching ultrasonic pulses and sampling the receivedwaveform as illustrated in FIG. 5 is repeated one or more times for eachelement in the array before proceeding to the next element. Then thesampled waveforms thus produced can be digitally averaged to produce acomposite waveform that will have a lesser random noise component thanany single sampled waveform. The number of repetitions necessary tosufficiently reduce noise can be determined by testing in a fashionknown to one skilled in the art.

Having thus reviewed the internal operation of the densitometer of FIGS.1-4, it is now possible to understand the methods of use of thedensitometer to measure the physical properties of the member. The firstmethod of use involves measuring transit time of an ultrasonic pulsethrough a subject and comparing that time to the time an ultrasonicpulse requires to travel an equal distance in a substance of knownacoustic properties such as water. To use the densitometer in thisprocedure, the adjustable arm 16 is adjusted until the member of thesubject, such as the heel, is clamped between the transducers 21. Thenthe knob 19 is tightened to fix the adjustable arm in place. Theactuator button 12 is then pressed to initiate a pulse and measurement.Next the densitometer is removed from the subject while keeping the knob19 tight so that the distance between the transducers 21 remains thesame. The device 10 is then placed about or immersed in a standardmaterial 32 with known acoustic properties, such as by immersion in abath of distilled water. The actuator button 12 is pressed again so thatacoustic signals are transmitted from the-transmit transducer 21 throughthe material 32 to the receive transducer 21. While it is advantageousto utilize the whole array of elements a through l for the measurementof the member, it may only be necessary to use a single pair of elementsfor the measurement through the standard assuming only that the standardis homogeneous, unlike the member. The signal profiles received by thetwo measurements are then analyzed by the microprocessor 38 and thesignal processor 41. This analysis can be directed both to thecomparative time of transit of the pulse through the subject as comparedto the standard and to the characteristics of the waveform in frequencyresponse and attenuation through the subject as compared to thestandard.

Thus in this method the densitometer may determine the physicalproperties and integrity of the member 32 by both or either of two formsof analysis. The densitometer may compare the transit time of theacoustic signals through the member with the transmit time of theacoustic signals through the material of known acoustic properties,and/or the device 10 may compare the attenuation as a function offrequency of the broadband acoustic signals through the member 32 withthe attenuation of corresponding specific frequency components of theacoustic signals through the material of known acoustic properties. The“attenuation” of an acoustic signal through a substance is thediminution of the ultrasonic waveform from the propagation througheither the subject or the standard. The theory and experiments usingboth of these methods are presented and discussed in Rossman, P. J.,Measurements of Ultrasonic Velocity and Attenuation In The Human OsCalcis and Their Relationships to Photon Absorptiometry Bone MineralMeasurements (1987) (a thesis submitted in partial fulfillment of therequirements for the degree of Master of Science at the University ofWisconsin-Madison). Tests have indicated that there exists a linearrelationship between ultrasonic attenuation (measured in decibels) (dB))at specific frequencies, and those frequencies. The slope (dB/MHz) ofthe linear relationship, referred to as the broadband ultrasonicattenuation, is dependent upon the physical properties and integrity ofthe substance being tested. With a bone, the slope of the linearrelationship would be dependent upon the bone mineral density. Thusbroadband ultrasonic attenuation through a bone is a parameter directlyrelated to the quality of the cancellous bone matrix.

The microprocessor 38 may therefore be programmed so that the devicedetermines the physical properties and integrity of the member bycomparing either relative transit times and/or relative broadbandultrasonic attenuation through the member and a material of knownacoustic properties. When comparing the transit times, themicroprocessor 38 may be programmed most simply so that the electronics,having received the acoustic signals after they have been transmittedthrough the member, determines the “member” transit time of thoseacoustic signals through the member, and after the acoustic signals havebeen transmitted through the material of known acoustic properties,determines the “material” transit time of the acoustic signals throughthe material. These time periods may be measured most simply by countingthe number of clock pulses of known frequency emitted by the timer 43between the time of launching the pulse and the sensing of the receivedpulse at the A/D converter 42. The microprocessor 38 then makes amathematical “time” comparison of the member transit time to thematerial transit time and then relates that mathematical time comparisonto the physical properties and integrity of the member. The mathematicaltime comparison may be made by either determining a difference betweenthe member transit time and the material transit time, or by determininga ratio between the member transit time and the material transit time.

As a second method of using the densitometer, it may also determine thephysical properties and integrity of the member 32 by determining andcomparing the attenuation of the broadband frequency components of theacoustic signals through the member without reference to a materialhaving known acoustic properties. Using this method, the comparison ofvelocity to a standard is not necessary and absolute transit time of thepulse need not be calculated since it is attenuation that is measured.In such a mode, it is preferable that the transmit transducer 21transmits an acoustic signal which has a broad range of frequencycomponents, such as a simple ultrasonic pulse. In any case, the acousticsignal should have at least one specific frequency component.

In this attenuation comparison mode, the microprocessor 38 is programmedso that after the receive transducer 21 receives the acoustic signalstransmitted through the bone member 32, it determines the absoluteattenuation through the member 32 of the frequency component spectrum ofthe acoustic signals. It is to facilitate the measurement of attenuationthat the excitation amplifier circuit 55 and the receiver amplifier 59have amplification levels which may be digitally controlled. Bysuccessively varying the gain of the amplifiers 55 and 59 on successivepulses, the circuit of FIG. 4 can determine what level of gain isnecessary to place the peak of the received waveform at a proper voltagelevel. This gain is, of course, a function of the level of attenuationof the acoustic pulse during transit through the member 32. After thereceive transducer 21 receives acoustic signals, microprocessor 38 inconjunction with the signal processor 41 determines the absoluteattenuation of individual specific frequency components of the receivedacoustic signal transmitted through the material. The digital signalprocessor 41 then makes mathematical “attenuation” comparisons of thecorresponding individual specific frequency components through themember. A set of mathematical attenuation comparisons betweencorresponding frequency components may be thereby obtained, onecomparison for each frequency component compared. The manner in whichthe attenuation functions with respect to frequency can thus be derived.The microprocessor 38 and digital signal processor 41 then relate thatfunction to the physical properties and integrity of the member.

Shown in FIG. 7 is a sample broadband ultrasonic pulse and a typicalreceived waveform. To achieve an ultrasonic signal that is very broad inthe frequency domain, i.e., a broadband transmitted signal, anelectronic pulse such as indicated at 70 is applied to the selectedultrasonic transducer in the transmit array 21 which then resonates witha broadband ultrasonic emission. The received signal, such as indicatedat 72 in FIG. 7 in a time domain signal plot, is then processed bydiscrete Fourier transform analysis so that it is converted to thefrequency domain. Shown in FIG. 8 is a pair of plots of sample receivedsignals, in frequency domain plots, showing the shift in received signalintensity as a function of frequency between a reference object and aplug of neoprene placed in the instrument. FIG. 9 illustrates a similarcomparison, with FIG. 8 using relative attenuation in the verticaldimension and FIG. 9 using power of the received signal using a similarreference material. Both representations illustrate the difference inrelative intensities as a function of frequency illustrating howbroadband ultrasonic attenuation varies from object to object. Theactual value calculated, broadband ultrasonic attenuation, is calculatedby first comparing the received signal against the reference signal,then performing the discrete Fourier transform to convert to frequencydomain, then performing a linear regression of the difference inattenuation slope to derive broadband ultrasonic attenuation.

The mathematics of the discrete Fourier transform are such that anotherparameter related to bone member density may be calculated in additionto, or in substitution for, broadband attenuation (sometimes referred toas “attenuation” or “BUA” below). When the discrete Fourier transform isperformed on the time-domain signal, the solution for each pointincludes a real member component and an imaginary member component. Thevalues graphed in FIGS. 8 and 9 are the amplitude of the received pulseas determined from this discrete Fourier transform by taking the squareroot of the sum of the squares of the real component and the imaginarycomponent. The phase angle of the change in phase of the ultrasonicpulse as it passed through the member can be calculated by taking thearctangent of the ratio of the imaginary to the real components. Thisphase angle value is also calculated to bone member density.

The microprocessor 38 may also be programmed so that the densitometersimultaneously performs both functions, i.e. determines both transittime and absolute attenuation of the transmitted acoustic signals, firstthrough the member and then through the material with known acousticproperties. The densitometer may then both derive the broadbandultrasonic attenuation function and make a mathematical time comparisonof the member transit time to the material transit time. Themicroprocessor 38 and digital signal processor 41 then relate both thetime comparison along with the attenuation function to the physicalproperties and integrity, or density of the member 32.

In yet another possible mode of operation, the microprocessor 38 may beprogrammed so that the densitometer 10 operates in a mode whereby theneed for calculating either the relative transit time or the attenuationof the acoustic signals through a material of known acoustic propertiesis eliminated. In order to operate in such a mode, the microprocessor 38would include a database of normal absolute transit times which arebased upon such factors as the age, height, weight, race or the sex ofthe individual being tested as well as the distance between thetransducers or the thickness or size of the member. This database ofnormal transit times can be stored in the non-volatile memory or couldbe stored in other media. When testing an individual in this mode, therelevant factors for the individual are placed into the microprocessor38 to select the pertinent normal transit time based on those factors.The transducers 21 are placed on the bone member being tested asdescribed above. When the actuator button 12 is pressed, the acousticsignals are transmitted through the member 32. The receive transducer 21receives those signals after they have been transmitted through themember, and the electronics 31 then determine the “member” transit timeof the acoustic signals through the member. The microprocessor 38 anddigital signal processor 41 then make a mathematical comparison of themeasured member transit time to the selected database normal transittime, and relate the mathematical time comparison to the physicalproperties and integrity, or density of the member, which is displayed.

As an alternative output of the densitometer of the present invention,the digital display 18 could also include a display corresponding to thepattern of the array of elements on the face of the transducer 21 asseen in FIG. 3. This display could then display, for each element athrough 1, a gray scale image proportional to the parameter, i.e.transit time or attenuation, being measured. This image may provide avisual indication to an experienced clinician as to the physicalproperties of the member present in the patient.

Shown in FIG. 6 is a circuit schematic for an alternative embodiment ofan ultrasonic densitometer constructed in accordance with the presentinvention. In the circuit of FIG. 6, parts having similar structure andfunction to their corresponding parts in FIG. 4 are indicated withsimilar reference numerals.

The embodiment of FIG. 6 is intended to function with only a singletransducer array 21 which functions both as the transmit and the receivetransducer array. An optional reflecting surface 64 may be placed on theopposite side of the member 32 from the transducer array 21. A digitallycontrolled multiple pole switch 66, preferably an electronic switchrather than a mechanical one, connects the input to and output from theelements of the transducer array 21 selectively either to the excitationamplifier 55 or to the controllable gain receiver/amplifier circuit 59.The switch 66 is connected by a switch control line 68 to an output ofthe microprocessor 38.

In the operation of the circuit of FIG. 6, it functions in most respectslike the circuit of FIG. 4, so only the differences need be discussed.During the launching of an ultrasonic pulse, the microprocessor 38causes a signal to appear on the switch control line 68 to cause theswitch 66 to connect the output of the excitation amplifier 55 to theselected element in the transducer array 21. Following completion of thelaunching of the pulse, the microprocessor 38 changes the signal on theswitch control line 68 to operate the switch 66 to connect the selectedelement or elements as an input to the amplifier 59. Meanwhile, thepulse propagates through the member 32. As the pulse transits throughthe member, reflective pulses will be generated as the pulse crossesinterfaces of differing materials in the member and, in particular, asthe pulse exits the member into the air at the opposite side of themember. If the transition from the member to air does not produce asufficient reflective pulse, the reflecting surface 64 can be placedagainst the opposite side of the member to provide an enhanced reflectedpulse.

The embodiment of FIG. 6 can thus be used to analyze the physicalproperties and integrity of a member using only one transducer 21. Allof the methods described above for such measurements may be used equallyeffectively with this version of the device. The transit time of thepulse through the member can be measured simply by measuring the timeperiod until receipt of the reflected pulse, and then simply dividing bytwo. This time period can be compared to the transit time, over asimilar distance, through a standard medium such as water. The timeperiod for receipt of the reflected pulse could also be simply comparedto standard values for age, sex, etc. Attenuation measurements to detectdifferential frequency measurement can be directly made on the reflectedpulse. If no reflecting surface 64 is used, and it is desired todetermine absolute transit time, the thickness of the member or samplecan be measured.

The use of the multi-element ultrasonic transducer array for thetransducers 21, as illustrated in FIG. 3, enables another advantageousfeature of the instrument of FIGS. 1-9. In using prior artdensitometers, it was often necessary to precisely position theinstrument relative to the body member of the patient being measured tohave useful results. The difficulty arises because of heterogeneities inthe bone mass and structure of actual body members. A measurement takenat one location of density may be significantly different from ameasurement taken close by. Therefore prior art instruments fixed thebody member precisely so that the measurement could be taken at theprecise location each time.

The use of the ultrasonic transducer array obviates the need for thisprecise positioning. Using the instrument of FIGS. 1-9, the instrumentperforms a pulse and response, performs the discrete Fourier transform,and generates a value for broadband ultrasonic attenuation for each pairof transducer elements a through l. Then the microprocessor 38 analyzesthe resulting array of bone ultrasonic density measurements toreproducibly identify the same region of interest each time. In otherwords, since the physical array of transducers is large enough toreliably cover at least the one common region of interest each time, themeasurement is localized at the same locus each time by electricallyselecting the proper location for the measurement from among thelocations measured by the array. The instrument of FIGS. 1-9 isconveniently used by measuring the density of the os calcis as measuredthrough the heel of a human patient. When used in this location, it hasbeen found that a region of interest in the os calcis can be locatedreliably and repeatedly based on the comparisons of broadband ultrasonicattenuation at the points in the array. The region of interest in the oscalcis is identified as a local or relative minimum in broadbandultrasonic attenuation and/or velocity closely adjacent the region ofhighest attenuation values in the body member. Thus repetitivemeasurements of the broadband ultrasonic attenuation value at this sameregion of interest can be reproducibly taken even though thedensitometer instrument 10 is only generally positioned at the samelocation for each successive measurement.

This technique of using a multiple element array to avoid positioncriticality is applicable to other techniques other than thedetermination of broadband ultrasonic attenuation as described here. Theconcept of using an array and comparing the array of results todetermine measurement locus would be equally applicable to measurementstaken of member-density based on speed of sound transit time, othermeasurements of attenuation or on the calculation of phase anglediscussed above. The use of such a multiple-element array, withautomated selection of one element in the region of interest, can alsobe applied to other measurement techniques useful for generatingparameters related to bone member density, such as measuring speedchanges in the transmitted pulse such as suggested in U.S. Pat. No.4,361,154 to Pratt, or measuring the frequency of a “sing-around”self-triggering pulse as suggested in U.S. Pat. No. 3,847,141 to Hoop.The concept which permits the position independence feature is that ofan array of measurements generating an array of data points from which aregion of interest is selected by a reproducible criterion or severalcriteria. The number of elements in the array also clearly can be variedwith a larger number of elements resulting in a greater accuracy inidentifying the same region of interest.

In this way, the ultrasound densitometer of the present inventionprovides a device capable of rapid and efficient determination of thephysical properties of a member in vivo without the use of radiation.Because the densitometer is constructed to operate under the control ofthe microprocessor 38, it can be programmed to operate in one of severalmodes, as discussed above. This allows both for flexibility to clinicalgoals as well as efficient use of the device.

Basin Embodiment

Shown in FIG. 10 is another variation on an ultrasonic densitometerconstructed in accordance with the present invention. In thedensitometer 100 of FIG. 10, there are two ultrasonic transducer arrays121, which are generally similar to the ultrasonic transducer arrays 21of the embodiment of FIG. 1, except that the transducer arrays 21 arefixed in position rather than movable.

The densitometer 100 includes a generally box-shaped mounting case 101with sloping upper face 102 in which is formed a basin 103. The basin103 is sized to receive a human foot and is generally trigonous along avertical plane aligned with the length of the foot so that when the footis placed within the basin 103, the toes of the foot are slightlyelevated with respect to the heel of the foot.

The transducer arrays 121 are positioned in the case 101 so that theyextend into the basin 103 to be on opposite sides of the heel of thefoot placed in the basin 103. When the foot is in position within thebasin 103, the sole of the foot may rest directly on a bottom 104 of thebasin 103 with the heel of the foot received within a curved pocket 106forming a back wall of the basin 103. As so positioned, the transducerarrays 121 are on either side of the os calcis. It has been demonstratedthat placing the transducer approximately 4 centimeters up from the soleand 3.5 centimeters forwardly from the rearward edge of the heel placesthe transducers in the desired region and focused on the os calcis.

The foot may, alternatively, rest on a generally planar foot plate 108having a contour conforming to the bottom 104 and placed against thebottom 104 between the foot and the bottom 104. The foot plate 108 holdsan upwardly extending toe peg 110 for use in reducing motion of the footduring the measurement process. Referring to FIG. 11, the toe peg 110 issized to fit between the big toe and the next adjacent toe of a typicalhuman foot and is mounted in a slot 112 so as to be adjustable generallyalong the length of the foot to accommodate the particular length of thefoot.

The slot 112 cants inward toward a medial axis 114 of the foot, definedalong the foot's length, as one moves along the slot 112 towards theportion of the foot plate 108 near the heel of the foot. This cantingreflects the general relation between foot length and width and allowssimple adjustment for both dimensions at once.

The toe peg 110 is sized to fit loosely between the toes of the footwithout discomfort and does not completely prevent voluntary movement ofthe foot. Nevertheless, it has been found that the tactile feedback tothe patient provided by the toe peg 110 significantly reduces footmovement during operation of the densitometer 100. Two different footplates 108, being mirror images of each other, are used for the left andright foot.

Referring to FIG. 12, the toe peg 110 is held to the slot 112 by afastener 111 having a threaded portion which engages correspondingthreads in the toe peg 110. The head of the threaded fastener 111engages the slot 112 so as to resist rotation. Thus, the toe peg 110 maybe fixed at any position along the length of the slot 112 by simplyturning the toe peg 110 slightly about its axis to tighten the threadedfastener 111 against the foot plate 108.

Referring again to FIG. 10, the basin 103 of the densitometer 110 isflanked, on the upper face 102 of the enclosure 101, by two foot restareas 116 and 118 on the left and right side respectively. Forexamination of a patient's right foot, the patient's left foot may reston foot rest area 118 while the patient's right foot may be placedwithin basin 103. Conversely, for examination of the patient's leftfoot, the left foot of the patient is placed within basin 103 and thepatient's right foot may rest on foot rest area 116. The foot rest areashave a slope conforming to that of the upper face 102 and approximatelythat of bottom 104. The flanking foot rest areas 116 and 118 allow thedensitometer 100 to be used in comfort by a seated patient.

When the densitometer 100 is not in use, the basin area 103 is coveredwith a generally planar cover 120 hinged along the lower edge of thebasin 103 to move between a closed position substantially within theplane of the upper face 102 and covering the basin 103, and an openposition with the plane of the cover 120 forming an angle a with thebottom 104 of the basin 103 as held by hinge stops 122. The angle a isapproximately 90° and selected so as to comfortably support the calf ofthe patient when the patient's foot is in place within basin 103. Tothat end, the upper surface of the cover 120, when the cover 120 is inthe open position, forms a curved trough to receive a typical calf.

The support of the patient's calf provided by the cover 120 has beenfound to reduce foot motion during operation of the densitometer 100.

Referring now to FIGS. 10 and 12, because the densitometer 100 employsfixed transducers 121, a coupling liquid is provided in the basin 103 toprovide a low loss path for acoustic energy between the transducers 121and the patient's foot regardless of the dimensions of the latter. Thecoupling liquid is preferably water plus a surfactant, the latter whichhas been found to improve the signal quality and consistency of thereading of the densitometer. The surfactant may be, for example, acommercially available detergent. It will be recognized, however, thatother flowable, acoustically conductive media may be used to provideacoustic coupling, and hence, that the term “coupling liquid” should beconsidered to embrace materials having a viscosity higher than that ofwater such as, for example, water based slurries and thixotropic gels.

For reasons of hygiene, the exhaustion of the surfactant, and possiblereduction of signal quality with the collection of impurities in thecoupling liquid, it has been determined that the liquid in the basin 103should be changed in between each use of the densitometer 103. Changingthis liquid is time consuming and ordinarily would require convenientaccess to a sink or the like, access which is not always available.Failure to change the liquid may have no immediate visible effect, andhence changing the liquid is easy to forget or delay. For this reason,the present embodiment employs an automated liquid handling systemlinked to the ultrasonic measurement operation through circuitrycontrolled by microprocessor 38 to be described.

Referring to FIG. 13 in the present embodiment, premixed water andsurfactant for filling the basin 103 are contained in a removablepolypropylene supply tank 124, whereas exhausted water and surfactantfrom the basin 103 are received by a similar drain tank 126. Each tank124 and 126 contains a manual valve 128 which is opened when the tanksare installed in the densitometer 100 and closed for transporting thetanks to a remote water supply or drain. The supply tank 124 and thedrain tank 126 have vents 150, at their upper edges as they are normallypositioned, to allow air to be drawn into or expelled from the interiorof the tanks 124 and 126 when they are in their normal position withinthe densitometer 100 and valves 128 are open. The tanks 124 and 126 holdsufficient water for approximately a day's use of the densitometer 100and thus eliminate the need for convenient access to plumbing.

The valve 128 of the supply tank 124 connects the tank through flexibletubing to a pump 130 which may pump liquid from the supply tank 124 to aheating chamber 132.

Referring to FIG. 14, the heating chamber 132 incorporates a resistiveheating element 164 which is supplied with electrical current through athermal protection module in thermal contact with the coupling liquid inthe heating chamber 132. The thermal protection module 166 includes athermostat and a thermal fuse, as will be described below. A thermistor168, also in thermal communication with the liquid in the heatingchamber, provides a measure of the liquid's temperature during operationof the densitometer 100. The heater chamber 132 additionallyincorporates an optical level sensor 172. The level sensor 172 detectsthe level of liquid in the heating chamber 132 by monitoring changes inthe optical properties of a prism system when the prism is immersed inliquid as opposed to being surrounded by air. The operation of thethermistor 168 and the level sensor 172 will be described further below.

Referring again to FIG. 13, the heating chamber 132 communicates throughan overflow port 134 and flexible tubing to an overflow drain outlet136. The overflow outlet 136 is positioned at the bottom of thedensitometer 100 removed from its internal electronics. The overflowport 134 is positioned above the normal fill height of the heatingchamber 132 as will be described in detail below.

The heating chamber 132 also communicates, through its lowermost point,with an electrically actuated fill valve 138 which provides a path,through flexible tubing, to a fill port 140 positioned in the wall ofbasin 103.

In the opposite wall of the basin 103 is an overflow port 142 whichopens into the basin 103 at a point above the normal fill height of thebasin 103 and which further communicates, through a T-connector 144, tothe drain tank 126.

A drain 146, in the bottom 104 of the basin 103, provides a path to anelectronically actuated drain valve 148. The drain valve 148 operates toallow liquid in the basin 103 to flow through the drain 146 to theT-connector 144 and into the drain tank 126. The overflow port 142 anddrain 146 incorporate screens 152 to prevent debris from clogging thetubing or the drain valve 148 communicating with the drain tank 126.

Referring now to FIGS. 10 and 13, the supply tank 124 and the drain tank126 are positioned within the case 101 of the densitometer 100 andlocated at a height with respect to the basin 103 so that liquid willdrain from the basin 103 into the drain tank 126 solely under theinfluence of gravity and so that gravity alone is not sufficient to fillthe basin 103 from supply tank 124 when fill valve 138 is open. Further,the heating chamber 132 is positioned above the basin 103 so that oncethe heating chamber 132 is filled with liquid by pump 130, the fillingof the basin 103 from the heating chamber 132 may be done solely by theinfluence of gravity. Accordingly, the operation of the densitometer infilling and emptying the basin 103 is simple and extremely quiet.

In those situations where plumbing is readily accessible, either or bothof the supply and drain tanks 124 and 126 may be bypassed and directconnections made to existing drains or supply lines. Specifically, thepump 130 may be replaced with a valve (not shown) connecting the heatingchamber 132 to the water supply line. Conversely, the connection betweenthe T-connector 144 and the drain tank 126 may re-routed to connect theT-connector 144 directly to a drain.

Even with the constant refreshing of the coupling liquid in the basin103 by the liquid handling system of the present invention, the liquidcontacting surfaces of the basin 103, the heating chamber 132, thevalves 138 and 148, and the connecting tubing are susceptible tobacterial colonization and to encrustation by minerals. The coatings ofcolonization or encrustation are potentially unhygienic andunattractive. Sufficient build-up of minerals or bacteria may alsoadversely affect the operation of the densitometer 100 either byrestricting liquid flow through the tubing, by interfering with theoperation of the valves 138 or 148, or by adversely affecting theacoustical properties of the transducer array 121.

For this reason, the densitometer 100 is desirably periodically flushedwith an antibacterial solution and a weak acid, the latter to removemineral build-up. These measures are not always effective or may beforgotten, and hence, in the present invention critical water contactingsurfaces are treated with a superficial antibacterial material which isalso resistant to mineral encrustation. The preferred treatment is theSPI-ARGENTÔ surface treatment offered by the Spire Corporation ofBedford, Mass. which consists of an ion beam assisted deposition ofsilver into the treated surfaces. The resulting thin film isbactericidal, fungistatic, biocompatible, and mineral resistant. Theproperties of being both bactericidal and fungistatic are generallytermed infection resistant.

This surface treatment is applied to the water contacting surfaces ofthe basin 103, the heating chamber 132 and the critical movingcomponents of the valves 138 and 148.

Referring now to FIG. 14, the general arrangement of the electricalcomponents of FIG. 4 is unchanged in the ultrasonic densitometer 100 ofFIG. 10 except for the addition of I/O circuitry and circuitry tocontrol the pump 130, valves 138 and 148, and heating chamber 132 of theliquid handling system. In particular, microprocessor 38 nowcommunicates through bus 40 with an I/O module 174, a pump/valve controlcircuit 160 and a heater control circuit 162.

I/O module 174 provides the ability to connect a standard video displayterminal or personal computer to the densitometer 100 for display ofinformation to the user or for subsequent post processing of the dataacquired by the densitometer and thus allows an alternative tomicroprocessor 38 and display 18 for processing and displaying theacquired ultrasound propagation data.

The pump/valve control circuit 160 provides electrical signals to thefill valve 138 and the drain valve 148 for opening or closing each valveunder the control of the microprocessor 38. The pump/valve controlcircuit 160 also provides an electrical signal to the pump 130 to causethe pump to begin pumping water and surfactant from the supply tank 124under the control of microprocessor 38, and receives the signal from thelevel sensor 172 in the heating chamber 132 to aid in the control of thepump 130 and valve 138.

The heater control circuit 162 controls the current received by theresistive heating element 164 and also receives the signal from athermistor 168 in thermal contact with the heating chamber 132. A secondthermistor 170, positioned in basin 103 to be thermal contact with theliquid in that basin 103, is also received by the heater control circuit162.

Referring now to FIGS. 13 and 14, during operation of the densitometer100 and prior to the first patient, the basin 103 will be empty, thesupply tank 124 will be filled and contain a known volume of water andsurfactant, and the drain tanks 126 will be empty. Both manual valves128 will be open to allow flow into or out of the respective tanks 124and 126 and the electrically actuated fill valve 138 and drain valve 148will be closed.

Under control of microprocessor 38, the pump/valve control circuit 160provides current to the pump 130 which pumps water and surfactant upwardinto heating chamber 132 until a signal is received from level sensor172. When the heating chamber 132 is filled to the proper level asindicated by level sensor 172, the signal from level sensor 172 topump/valve control circuit 160 causes the pump 130 to be turned off. Atthis time, a predetermined volume of liquid is contained in heatingchamber 132 which translates to the proper volume needed to fill basin103 for measurement.

Under command of microprocessor 38, the heater control circuit 162provides a current through thermal protection module 166 to resistiveheating element 164. The temperature of the liquid in the heatingchamber 132 is monitored by thermistor 168 and heating continues untilthe liquid is brought to a temperature of approximately 39° C. Thethermistor and a thermal fuse (not shown) of the thermal protectionmodule 166 provide additional protection against overheating of theliquid. The thermistor opens at 50° C. and resets automatically as itcools and the thermal fuse opens at 66° C. but does not reset and mustbe replaced. The opening of either the thermistor or the thermal fuseinterrupts current to the resistive heating element 164.

When the liquid in the heating chamber 132 is brought to the correcttemperature, fill valve 138 is opened by microprocessor 38, throughpump/valve control circuit 160, and liquid flows under the influence ofgravity into the basin 103 at the proper temperature. The control of thetemperature of the liquid serves to insure the comfort of the patientwhose foot may be in the basin 103 and to decrease any temperatureeffects on the sound transmission of the water and surfactant.

Once the heated liquid has been transferred from the heating chamber 132to the basin 103, the fill valve 138 is closed and the pump 130 isreactivated to refill the heating chamber 132. Thus, fresh liquid forthe next measurement may be heated during the present measurement toeliminate any waiting between subsequent measurements.

With liquid in place within the basin 103, the measurement of the oscalcis by the densitometer 100 may begin. In this respect, the operationof the ultrasonic densitometer of FIG. 10 is similar to that of theembodiment of FIG. 1 except that the order of pulsing and measurementcan be varied. In the apparatus of FIG. 1, the measurement pulse throughthe member was generally performed before the reference pulse throughhomogenous standard, i.e. water. In the densitometer 100 of FIG. 10,since the distance between the transducers 121 is fixed, the referencepulse through the homogenous standard material, which is simply theliquid in basin 103, may be conducted before or after a measurementpulse through a live member is performed. In fact, because thetemperature of the liquid in the basin 103 is held steady by thetemperature control mechanism as described, the standard transmit timemeasurement can be made once for the instrument and thereafter onlymeasurement pulses need be transmitted.

Preferably, the standard transit time measurement is stored as a numberin the memory of microprocessor 38 during the initial calibration of theunit at the place of manufacture or during subsequent recalibrations.During the calibration of the densitometer 100, the signal from thethermistor 170 is used to produce a transit time corrected for thetemperature of the liquid according to well known functional relationslinking the speed of sound in water to water temperature. It is thiscorrected transit time that is stored in the memory associated withmicroprocessor 38 as a stored standard reference.

The transit time of the measurement pulses is compared to the storedstandard reference transit times through the coupling liquid to give anindication of the integrity of the member just measured. Thus, one maydispense with the reference pulse entirely. Empirical tests havedetermined that by proper selection of a standard reference value storedin the memory of microprocessor 38 and by holding the liquid in thebasin within a temperature range as provided by the heating chamber 132,no reference pulse need be launched or measured.

Using this variation, a mathematical comparison of the measured transittime, or transit velocity, must be made to the standard. Since, in theinterests of accuracy, it is preferred to use both changes in transittime (velocity) and changes in attenuation to evaluate a member in vivo,the following formula has been developed to provide a numerical valueindicative of the integrity and mineral density of a bone:

bone integrity value=A(SOS-B)+C(BUA-D)  (1)

In this formula, “SOS” indicates the speed of sound, or velocity, of themeasurement ultrasonic pulse through the member, and is expressed inmeters per second. The speed of sound (SOS) value is calculated from themeasured transit time by dividing a standard value for the member widthby the actual transit time measured. For an adult human heel, it hasbeen found that assuming a standard human heel width of 40 mm at thepoint of measurement results in such sufficient and reproducibleaccuracy that actual measurement of the actual individual heel is notneeded.

BUA is broadband ultrasonic attenuation, as described in greater detailabove. The constants A, B, C, and D offset and scale the influence ofthe BUA measurement relative to the SOS measurement to provide a moreeffective predictor of bone density. These constants may be determinedempirically and may be selected for the particular machine to providenumbers compatible with dual photon absorptiometry devices and to reducebone width effects. Since this method utilizing ultrasonic measurementof the heel is quick and free from radiation, it offers a promisingalternative for evaluation of bone integrity.

The densitometer 100 may be used with or without an array of ultrasonicransducers in the transducers 121. In its simplest form the mechanicalalignment of the heel in the device can be provided by the shape andsize of the basin 103. While the use of an array, and region-of-interestscanning as described above, is most helpful in ensuring a reproducibleand accurate measurement, mechanical placement may be acceptable forclinical utility, in which case only single transducer elements arerequired.

Upon completion of the measurement, the drain valve 148 is opened bymicroprocessor 38 through pump/valve control circuitry 160, and theliquid in the basin 103 is drained through “T” 144 to the drain tank126. At the beginning of the next measurement, the drain valve 148 isclosed and liquid is again transferred from the heating chamber 132 ashas been described.

With repeated fillings and drainings of the basin 103, the level ofliquid in the fill tank 124 decreases with a corresponding increase inthe level of the liquid in the drain tank 126. The height of the liquidin each tank 124 and 126 may be tracked by a conventional level sensorsuch as a mechanical float or a capacitive type level sensor.

Preferably no additional level sensor is employed. The volume of liquidfor each use of the densitometer 100 is known and defined by the filllevel of the heating chamber 132. The microprocessor 38 may thereforetrack the level of liquid remaining in the supply tank 124 by countingthe number of times the basin 103 is filled to provide a signal to theuser, via the display 18 or a remote video display terminal (not shown),indicating that the tanks 124 and 125 need to be refilled and drainedrespectively. This signal to the user is based on the number of timesthe basin 103 is filled and a calculation of the relative volumes of theheating chamber 132 and supply tank 124.

After completion of the use of the densitometer 100 for a period oftime, the densitometer may be stored. In a storage mode, after both thesupply tank 124 and drain tank 126 have been manually emptied, themicroprocessor 38 instructs the pump/valve control circuit 160 to openboth the fill valve 138 and the drain valve 148 and to run the pump 130.The drain valve 138 is opened slightly before the pump 130 is actuatedto prevent the rush of air from causing liquid to flow out of theoverflow port 134.

Referring now to FIGS. 10 and 15, the transducers 121 are inserted intothe basin 103 through tubular sleeves 180 extending outward from thewalls of the basin 103 at the curved pocket along an axes 212 of theopposed transducers 121. The tubular sleeves 180 define a circular borein which the transducers 121 may be positioned. Each transducer 121seals the sleeve 180 by compression of o-ring 182 positioned on theinner surface of the sleeve 180.

Although the transducers 121 fit tightly within the sleeves 180, theirseparation and alignment are determined not by the sleeves 180 but by anindependent C-brace 184 comprising a first and second opposed arm 186separated by a shank 188. A transducers 121 is attached to one end ofeach of the arms 186, the other ends of the arms 186 fitting against theshank 188.

The arms 186 are generally rectangular blocks transversely bored toreceive the cylindrically shaped transducers 121 at one end and to holdthem along axis 212. The other ends of the arms 186 provide planar facesfor abutting the opposite ends of the block like shank 188, the abuttingserving to hold the arms 186 opposed and parallel to each other.

Although the angles of the arms 186 with respect to the shank 188 aredetermined by the abutment of the planar faces of the arms 186 and theends of the shank 188, alignment of the arms 186 with respect to theshank 188 is provided by dowel tubes 190 extending outward from each endof the shank 188 to fit tightly within corresponding bores in the firstand second arm 186.

Cap screws 194 received in counterbored holes in the arms 186 passthrough the arms 186, the dowel tubes 190 are received by threaded holesin the shank 188 to hold the arm 186 firmly attached to the shank 188.The dowel tubes 190 and surfaces between the arms 186 and shank 188serve to provide extremely precise alignment and angulation of thetransducers 121, and yet a joint that may be separated to permit removalof the transducers 121 from the densitometer 10 for replacement orrepair.

Transducers 121 are matched and fitted to the arms 186 in a controlledfactory environment to provide the necessary acoustic signal strengthand reception. In the field, the shank 188 may be separated from one orboth arms 186 by loosening of the cap screws 194 so as to allow thetransducers 121 extending inward from the arms 186 to be fit within thesleeves 180. Proper alignment and angulation of the transducers is thenassured by reattaching the arm or arms 186 removed from the shank 188 tothe shank 188 to be tightened thereto by the cap screws 194. Thus, thealignment of the transducers is not dependent on the alignment of thesleeves 180 which may be molded of plastic and thus be of relatively lowprecision. Nor must alignment be tested while the transducers are in thesleeves 180 attached to the basin 103 but may be checked in a centralcontrolled environment.

Flexible Bladder Embodiment

Referring now to FIGS. 16 and 17, in yet another embodiment of thepresent invention, the opposed transducers 121 are fitted with annularcollars 200 which in turn are attached to flexible bladders 202extending inward to the basin 103, each bladder 202 containing a liquidor semi-liquid coupling “gel” 204.

The bladders 202 serve to contain the gel about the face of thetransducers 121 and conform to the left and right sides of a patient'sheel 207, respectively, to provide a path between the transducers 121and the soft tissue and bone of the heel 207 without intervening air.The bladder 202 further prevents the coupling material from directcontact with the heel to permit selection of the coupling gel 204 from abroader range of materials.

Compression of the bladders 202 against the heel 207, so as to providethe necessary coupling, is provided by a telescoping shank 181 shown inFIG. 16. In this alternative embodiment of the C-brace 184 of FIG. 15,the shank 188′ has been cut into two portions 206 and 208 slidablyconnected together by dowel pins 210 to provide necessary motion of thetransducers 121 inward along their axis to compress the bladders 202against the heel 207. One end of each dowel pin 210 is press fit withinbores in the 20 shank 188′ parallel to the axis 212 of the opposedtransducers in portion 206. The other ends of the dowel pins 210 slidewithin larger bores in portion 208 so that portions 208 and 206 mayslide toward and away from each other parallel to the axis 212. Withsuch motion, the attached arms 186 move towards and away from each otheradjusting the separation of the transducers 121 between an open positionfor insertion of the heel 207 and a closed position of known separationand orientation where portions 208 and 206 abut.

Control of the separation is provided by means of cam pins 214protruding from portions 206 and 208 on the side away from the extensionof the arms 186 and generally perpendicular to the axis 212. These pins214 are received by spiral shaped slots in a cam disk 217 fitting overthe cam pins 214. The disk includes radially extending lever 218 whosemotion rotates the disk causing the cam pins 214 within the slots 215 tobe moved together or apart depending on motion of lever 218.

Thus, the transducers 121 may be moved apart together with the bladders202 for insertion of the heel 207 into the basin 103. Once the heel isin place, motion of the lever 218 closes the transducers 121 to apredetermined fixed separation compressing the bladders 202 snuglyagainst the sides of the heel 207. The elasticity of the bladder filledwith coupling gel 204 provides an expanding force against the heel 207to closely conform the surface of the bladder 202 to the heel 207.

Cancellation of Heel Width Variations

Referring to FIGS. 17 and 18, generally the thicker the calcaneus 216 ofthe heel 207, the greater the attenuation of an acoustic signal passingthrough the heel 207 between transducers 121. Correspondingly, withgreater attenuation, the slope of attenuation as a function offrequency, generally termed broadband ultrasonic attenuation (BUA) alsoincreases as shown generally in FIG. 18 by plot 209. This assumesgenerally that the coupling medium 204 is of low or essentially constantattenuation as a function of frequency. Greater BUA is generallycorrelated to higher bone quality.

For constant heel thickness, lower TOF (faster sound speed) correspondsgenerally to higher bone quality. The time of flight (TOF) of anacoustic pulse between the transducers 121 will be proportional to thetime of flight of the acoustic pulse through regions A of FIG. 17comprising the path length through coupling gel 204, regions Bcomprising the path length through soft tissue of the heel 207surrounding the calcaneus 216, and region C comprising the path lengththrough the heel bone or calcaneus 216. Thus, $\begin{matrix}{{TOF} = {{\frac{1}{V_{A}}A} + {\frac{1}{V_{B}}B} + {\frac{1}{V_{C}}C}}} & (2)\end{matrix}$

where V_(A), V_(B), and V_(C) are the average speed of sound through thecoupling gel, soft tissue and bone respectively and A, B, C are the pathlengths through these same materials. Provided that the separationbetween the transducers 121 is a constant value K, then time of flightwill equal: $\begin{matrix}{{TOF} = {{\frac{1}{V_{A}}\left( {K - C - B} \right)} + {\frac{1}{V_{B}}B} + {\frac{1}{V_{C}}C}}} & (3)\end{matrix}$

The change in time of flight as a function the thickness of the bone C(the derivative of TOF with respect to C) will thus generally be equalto: $\frac{1}{V_{C}} - {\frac{1}{V_{A}}.}$

Referring now to FIG. 18, if the velocity of sound through the couplingmedium 204 is greater than that through the bone being measured$\left( {{V_{A} > V_{C}},\quad {{{or}\quad \frac{1}{V_{C}}} > \frac{1}{V_{A}}}} \right),$

then the functional relationship of TOF to heel width will be one ofincreasing as the heel becomes wider (indicated at plot 213 showingvalues of 1/TOF). On the other hand, if the velocity of sound throughthe coupling medium 204 is less than that through the bone beingmeasured$\left( {{V_{C} > V_{A}},{{{but}\quad \frac{1}{V_{A}}} > \frac{1}{V_{C}}}} \right),$

then the functional relationship of TOF to heel width will be one ofdecreasing as the heel becomes wider (indicated at plot 211 showingvalues of 1/TOF).

A combined bone health figure may be obtained by combining BUA and 1/TOFmeasurements (1/TOF because BUA increases but TOF decreases withhealthier bone). Further, if (1) the conditions of ultrasonicpropagation are adjusted so that the slope of 1/TOF with heel width isopposite in sign to the slope of BUA with heel width (i.e., V_(A)>V_(C))and (2) the BUA and 1/TOF measurements are weighted with respect to eachother so that the opposite slopes of the BUA and 1/TOF are equal, thenthe algebraic combination of the BUA and TOF, through addition forexample, will produce a bone quality measurement substantiallyindependent of heel width for a range of bone qualities.

This can be intuitively understood by noting that as the heel getswider, it displaces some of the coupling gel 204 from between the heel207 and each transducer 121, and by displacing material that conductssound slower than the bone being measured increasing the total speedwith which the sound is conducted.

Note that a similar effect may be obtained by proper scaling andcombination of BUA and TOF by multiplication and that other functions ofattenuation and TOF could be used taking advantage of their functionalindependence and their functional dependence in part on heel width.

Referring now to FIG. 19, generally BUA and TOF are functionally relatedto both bone quality and bone width. It should be possible, therefore,to solve the equations governing these relationships for bone qualityalone and thus to eliminate the effect of the common variable of heelwidth. With such an approach, the variable of heel width is eliminatednot just for a portion but through the entire range of bone measurementprovided that the coupling medium is different from the bone beingmeasured so that there will be a width effect in both BUA and TOFmeasurements.

Approximations of the algebraic relationships describing the functionaldependence of BUA and TOF on bone quality and bone width, can beobtained through the construction of a set of bone phantoms of differentwidths and bone qualities when using a particular coupling gel.Generally, for each value of BUA or TOF the data will describe a curve222 linking that value with different combinations of bone quality andbone width. This data may be placed in a look-up table in the memory ofthe microprocessor of the densitometer as has been previously described.

After BUA and TOF values are determined, the data of the look-up table(comprising many bone quality and bone width pairs for each of thedetermined BUA and TOF values) are scanned to find a bone quality andwidth data pair for the BUA value matching a bone quality and width datapair for the TOF value. This is equivalent to finding the intersectionof the two curves 222 associated with the measured BUA and TOF values.The matching bone quality values of the data base will give a bonequality having little or no bone width influence. This value may bedisplayed to the clinician. It is noted that the previously describedtechnique of summing weighted values of BUA and 1/TOF is but aspecialized form of this process of algebraic solution.

Alternatively, a matching bone width value can be identified, being thewidth of the measured heel, and used to correct either of the BUA or TOFvalues for display to the clinician in circumstances where BUA or TOFvalues are preferred for diagnosis.

This ability to cancel out heel width effects will work only for bonequalities where the relationship between the coupling gel 204 and thecalcaneus 216 are such as to provide a functional dependence on heelwidth. Cancellation will not occur, for example, if the density of thecalcaneus 216 being measured is substantially equal to the sound speedof the coupling gel 204 and thus where displacement of the coupling gelby similar bone will have no net effect on time of flight. Thus thecoupling gel must be properly selected. In this case, materials havinghigher sound speed may be selected for the coupling material. Thedifference between the coupling gel and the bone being measured willinfluence the accuracy of the cancellation of heel width effects.

Moderating this desire to improve heel width effects is the importanceof keeping the coupling gel 204 close to the acoustic properties of thesoft tissue of the heel 207 both to prevent reflection by impedancemismatch and to prevent variations in the thickness of the soft tissuein regions B from adding additional uncertainty to the measurement. Thecoupling medium of water provides good matching to the soft tissue ofthe heel 207 and has a sound velocity very close to bone and someosteoporotic conditions. Weighting of the attenuation and propagationtime may be made for water.

Although the preferred embodiment of the invention contemplates displayof a bone quality value or corrected TOF or BUA values, it will berecognized that the same effect might be had by displaying uncorrectedBUA or TOF values on a chart and establishing a threshold for healthy orweak bone based on the corrections determined as above.

Ultrasonic Densitometer with Scannable Focus

Referring now to FIG. 20, a receiving transducer array 300, similar toarray 21 described with respect to FIG. 1, may be positioned adjacent tothe heel of a patient (not shown) to receive an ultrasonic wave 410along axis 304. The receiving transducer array 300 includes apiezoelectric sheet 302 of substantially square outline positionednormal to the transmission axis 304 and is divided into transducerelements 400 as will be described, each which receives a differentportion of the ultrasonic wave 410 after passage through the heel.

The piezoelectric sheet 302 may be constructed of polyvinylidenefluoride and has a front face 306 covered with a grid of interconnectedsquare electrodes 308 deposited on the front face by vacuummetallization. These square electrodes 308 are arranged at theinterstices of a rectangular grid to fall in rectilinear rows andcolumns. Referring also to FIG. 23, each square electrode 308 is spacedfrom its neighboring electrodes 308 by approximately its width. Thesesquare electrodes 308 are connected together by metallized traces (notshown) and to a common voltage reference by means of a lead 310.

In the manufacture of the piezoelectric sheet 302, the polyvinylidenesheet is polarized to create its piezoelectric properties by heating andcooling the sheet in the presence of a polarizing electrical fieldaccording to methods generally understood in the art. In the presentinvention, this polarizing field is applied only to the area under thesquare electrodes 308 so that only this material is piezoelectric andthe material between square electrodes 308 has reduced or nopiezoelectric properties. As will be understood below, this selectivepolarization of the piezoelectric sheet 302 provides improved spatialselectivity in distinguishing between acoustic signals received atdifferent areas on the piezoelectric sheet.

Referring now to FIG. 22, opposite each electrode 308 on the back sideof the piezoelectric sheet 302 furthest from the source of theultrasonic wave 410 is a second electrode 312 having substantially thesame dimensions as the square electrodes 308 and aligned withcorresponding square electrodes 308 along transmission axis 304.

Referring to FIGS. 20 and 21, a connector board 318 of areal dimensionsubstantially equal to the piezoelectric sheet 302 has, extending fromits front surface, a number of conductive pins 320 corresponding to thepads 316 in number and location. The pins 320 are stake-type terminalsmounted to an epoxy glass printed circuit board 322 of a type well knownto those of ordinary skill in the art. Each conductive pin 320 isconnected directly to a preamplifer and then by means of printed circuittraces to a multiplexer 325 to a reduced number of control and datalines 324 which may be connected to the microprocessor 38 of thedensitometer through an A to D converter 42 described previously withrespect to FIG. 1 and as is well understood in the art. The preamplifersallow grounding of those electrodes 312 not active during scanning toreduce cross-talk between electrodes 312.

As shown in FIG. 21, the pins 320 of the connector board 318 areelectrically connected to electrodes 312 on the back surface of thepiezoelectric sheet 302 by means of a strip of thin (0.0005″) Mylar 316having conductive fingers 314 on its surfaces. The conductive fingers314 on the front and rear surfaces of the Mylar strip 316 are inelectrical communication through a plated-through hole 313 in the Mylar316 connecting the fingers 314.

Each conductive pin 320 is attached to a conductive finger 314 at oneedge of the Mylar strip 316 at the rear of the Mylar strip 316(according to the direction of the acoustic wave) by means of ananisotropically conductive adhesive film 315 providing electricalconduction only along its thinnest dimension, thus from pin 320 tofinger 314 but not between fingers 314 or pins 320. Anisotropicallyconductive film suitable for this purpose is commercially available from3M corporation of Minnesota under the trade name of 3M Z-Axis AdhesiveFilm.

The other end of each plated finger 314 on the front of the Mylar strip316 is then connected to an electrode 312 by a second layer ofanisotropically conductive adhesive film 317. The Mylar strip 316 flexesto allows the pins 320 to be spaced away from the electrode 312 toreduce reflections off the pins 320 such as may cause spurious signalsat the piezoelectric sheet 302. The Mylar strip 316 and conductivefingers 314 are essentially transparent to the acoustic wave.

Referring to FIG. 22, the Mylar strips 316 and adhesive film 315 and 317allow rapid assembly of the transducer 300. A single layer of conductivefilm 317 (not shown in FIG. 22) may be applied over the entire rearsurface of the piezoelectric sheet 302 and electrodes 312. Next aplurality of overlapping Mylar strips 316 may be laid down upon thissurface, each Mylar strip 316 extending laterally across thepiezoelectric sheet 302 with transversely extending conductive fingers314 for each electrode 312 of one row of conductive electrode 312. Theoverlapping of the Mylar strips 316 ensures that only a front edge ofeach strip 316 adheres to the piezoelectric sheet 302. Guide holes 319in the laterally extreme edges of the Mylar strips 316 fit into pins ina jig (not shown) to ensure alignment of fingers 314 with electrodes312.

Next, a second layer of the anisotropically conductive adhesive film 315is placed on the rear surfaces of the overlapping Mylar strips 316 andthe conductive pins 320 pressed down on this film 315, aligned with theother ends of the conductive fingers 314 to attach to their respectivefingers 314. The conductive pins 320 are then raised and fixed in spacedapart relationship with the piezoelectric sheet 302, the Mylar strips316 flexing to accommodate this displacement.

The ultrasonic wave 410 passing through portions of the piezoelectricsheet 302 between electrodes 308 and 312 may thereby be measured at anumber of points over the surface of the piezoelectric sheet by theelectric signals generated and collected by electrodes 308 and 312according to multiplexing methods well known in the art. Each electrodepair 308 and 312 provides an independent signal of the acoustic energypassing through the area of the piezoelectric sheet 302 embraced by theelectrode pair.

A protective frame 325 encloses the piezoelectric sheet 302 andconnector board 318 protecting them from direct contact with water ofthe basin 103 shown in FIGS. 10 and 15 into which the receivingtransducer array 300 may be placed. The frame 325 holds on its frontface an acoustically transparent and flexible material 326 such as aTeflon film so that the ultrasonic wave 410 may pass into the frame toreach the piezoelectric sheet 302.

The above described array may be used either to receive or transmitacoustic waves and is not limited to use in the medical area but mayprovide an inexpensive and rugged industrial acoustic array useful for avariety of purposes including industrial ultrasonic imaging and theconstruction of high frequency synthetic aperture microphones.

Positioned behind the frame 325 is an electric motor 328 driving acentral gear 330 about an axis aligned with transmission axis 304 andapproximately centered within the frame 325. The central gear 330 inturn engages two diagonally opposed planet gears 332 also turning aboutaxes aligned with the transmission axis. Each planet gear 332 has a rod334 extending forwardly from a front face of the planet gear 332 butoffset from the planet gear's axis to move in an orbit 336 thereabout.The orbit 336 has a diameter approximately equal to the spacing betweenelectrodes 308.

The rods 334 engage corresponding sockets 338 on the back side of theframe 325 at its opposed corners. Thus activation of the motor 328causes the piezoelectric sheet 302 and connector board 318 to follow theorbit 336 while maintaining the rows and columns of detector elements400 in horizontal and vertical alignment, respectively.

Referring now to FIG. 23, a sampling of the signals from the detectorelements 400 may be made at four points 342 in the orbit 336 at whicheach electrode 308 is first at a starting position, and then is movedhalf the inter-electrode spacing upward, leftward, or upward andleftward. The effect of this motion of the detector elements 400 is todouble the spatial resolution of the received acoustic signals withoutincreasing the amount of wiring or the number of detector elements 400.The sampling of acoustic energy at each of the points 342 is stored inthe memory of the microprocessor and can be independently processed toderive attenuation, BUA or time of flight measurements or a combinationof these measurements. These measurements are then converted to anintensity value of an image so that each pixel of the image has anintensity value proportional to the measured parameter. A clinicianviewing the image thus obtains not merely an image of the bone, but animage that indicates bone quality at its various points.

A transmitting ultrasonic transducer 408 is positioned opposite thereceiving transducer array 300 from the heel 207 and produces agenerally planar ultrasonic wave 410 passing into the heel. Generally,the acoustic signal received by each transducer element 400 will havearrived from many points of the heel.

Referring now to FIG. 24, if the transducer elements 400 were focused asindicated by depicted transducer elements 400′ to follow a hemisphere402 having a radius and hence focus at a particular volume element orvoxel 404 within the heel, acoustic signals from other voxels could becanceled providing greater selectivity in the measurement. In thisfocusing of the transducer elements 400′, the signals from each of theelements 400′ are summed together. Constructive and destructiveinterference of ultrasonic waves 410 from the heel 207 serve toeliminate acoustic signals not flowing directly from focus volumeelement 404.

For example as depicted, two acoustic signals 405 and 406 from focusvoxel 404 both crest at the location of a transducer element 400′ as aresult of the equidistance of each transducer element 400′ from focusvoxel 404. When the signals from transducer elements 400′ are summed,the signal from focus voxel 404 will increase. In contrast, acousticsignals from other voxels not equidistant to transducer elements 400′will tend to cancel each other when summed and thus decrease.

The present invention does not curve the transducer elements 400 into ahemisphere but accomplishes the same effect while retaining thetransducer element 400 in a planar array by delaying the signalsreceived by the transducer elements 400 as one moves toward thecentermost transducer element 400″ so as to produce an effectivehemispherical array. Like a hemispherical array, the center-mosttransducer elements 400 appear to receive the acoustic wave a littlelater than the transducer elements 400 at the edge of the receivingtransducer array 300. By using a phase delay of the signals instead ofcurving the receiving array 300, the position of the focus voxel 404 atwhich the receiving array 300 is focused, may be scanned electrically aswill be described. The signals from each of the transducer elements 400are received by the A/D converter 42 and stored in memory. Phaseshifting as described simply involves shifting the point at which onestarts reading the stored signals.

Adjusting the phase of the acoustic signals received by each of thetransducer elements 400 allows the location of the focus voxel 404 fromwhich data is obtained to be scanned through the heel. The phase issimply adjusted so that the effective arrival time of an acoustic signaloriginating at the desired location is the same for each of thetransducer elements 400.

Referring now to FIG. 25, the location of focus voxel 404 may be movedin a first and second raster scan pattern 412 and 414 (as readings aretaken over many ultrasonic pulses) to obtain separated planes of datanormal to the transmission axis 304. The first plane of data 412 may,for example, be positioned near the outer edge of the os calcis 216 tomeasure the cortical bone quality while the second plane 414 may beplaced in a centered position in the trabecular bone to obtain asomewhat different reading, both readings providing distinct data aboutthe bone.

It will be understood that this same approach of scanning in differentplanes may be used to obtain a volume of data within the heel 207, inthis case, the focus voxel 404 being moved to points on a threedimensional grid.

In another embodiment (not shown) the transmitting ultrasonic transducermay be an array and the phases of the ultrasonic signals transmitted byeach of the elements of the array may be phased so as to focus on aparticular voxel within the heel. In this case, the receiving array maybe a single broad area detector or may also be an array focused on thesame voxel for increased selectivity. The focus point of thetransmitting and receiving arrays may also be shifted with respect toeach other to investigate local sound transfer phenomenon. As before,the focal points of either array may be steered electrically by themicroprocessor through a shifting of the phases of the transmitted andreceived signals. To collect data, each element of the transmit arraymay be energized individually while all receive elements of the receivearray are read. This may be continued until each of the elements of thetransmit array have been energized.

Alternatively, referring to FIG. 28, the receiving array 300 may beactually formed so that its elements follow along the hemisphere 402 soas to have a fixed focus on focus voxel 404. Additional circuitry toeffect the phase adjustment needed to focus the array is not needed inthis case. The receiving array 300 is attached to an X-Y-Z table 600providing motion in each of three Cartesian axes under the control ofthe microprocessor via stepper motors 610. At each different location ofthe table 600, data may be collected from focus voxel 404 to establishthe data points on the three dimensional grid. The transmitting array408 may be held stationary or may be moved with the scanning of thereceiving array 300 and may be focused as well.

Referring now to FIG. 26, such a data volume 415 may include a pluralityof data voxels 416 each providing a measured member parameter for thebone or tissue at that point in the heel. A point of minimum bonedensity 418 may be found within this data volume 415 and used toidentify a region of interest 420 which will serve as a standard regionfor measuring the bone density of the heel. This region may beautomatically found after collection of the data volume 415 and onlythose voxels 416 within the region of interest 420 may be used for adisplayed measurement. This automatic location of a region of interest420 provides a much more precise bone characterization.

Acquiring a data volume 415 also provides the opportunity to use theextra data outside the region of interest 420 to ensure that the sameregion of interest 420 is measured in the patient's heel over a seriesof measurements made at different times. The data volume 415 may bestored in memory as a template that may be matched to subsequentlyacquired data volumes. The region of interest 420 spatially located withrespect to the first template, may then be used as the region ofinterest for the subsequent data volumes aligned with that template toprovide more repeatability in the measurement.

Referring now to FIG. 27 in such a template system in a first step 500,a collection of a data volume 415 within the heel is obtained. Atdecision block 502, if this is a first measurement of a particularpatient, a region of interest 420 is identified at process block 504from this data, as a predetermined volume centered about a point ofminimum bone density 418 as described with respect to FIG. 26. Atprocess block 506, the data volume is stored as a template along withthe region of interest defined with respect to the data of the template.

Referring again to decision block 502 on a subsequent measurement of apatient, the program may proceed to process block 508 and the templatepreviously established may be correlated to a new data volume 415collected at process block 500. The correlation process involvesshifting the relative locations of the two data volumes to minimize adifference between the values of each data voxel 416 of the datavolumes. In most situations, this will accurately align the two datavolumes so that corresponding voxels 416 of the two data volumes 415measure identical points within the patient's heel. The region ofinterest 420 associated with the template is then transferred to the newdata volume as it has been shifted into alignment with the template sothat the identical region of interest may be measured in a patient evenif the patient's foot has taken a different alignment with respect tothe transducer array 300 and 408. This use of the template's region ofinterest 420 is indicated by process block 510.

At process block 512, an index is calculated at the region of interest420 for the new data volume 415 being typically an average value of abone parameter such as BUA or time of flight for the voxels 416 withinthe region of interest 420. This data is then displayed to the clinicianat process block 520 as has been described.

It is specifically intended that the present invention not bespecifically limited to the embodiments and illustrations containedherein, but embrace all such modified forms thereof as come within thescope of the following claims.

We claim:
 1. A method of non-invasive and quantitative assessment of thestatus of bone tissue in vivo for at least one of the quantities,bone-mineral density, bone strength, bone fracture risk, bonearchitecture and bone quality comprising the steps of: positioning afirst and second transducer distant from the skin on opposite sides ofsaid bone tissue to provide gaps therebetween allowing free insertion ofsaid bone tissue unobstructed by the first and second transducers;acoustically coupling a first transducer and a second transducer to acompliant medium bridging the gaps between the first and secondtransducers and the distant skin; acoustically coupling the compliantmedium to the distant skin on opposite sides of said bone tissue whereinone of said first and second transducers is a single-element transducerand another of said first and second transducers is an array transducer;generating an ultrasound signal and directing said ultrasound signalfrom said first transducer to said second transducer through said bonetissue to obtain a bone-oriented output signal; and processing saidbone-oriented output signal whereby an estimate of said at least one ofthe quantities, bone-mineral density, bone strength, bone fracture risk,bone architecture and bone quality is obtained.
 2. The method of claim 1wherein said processing step includes determining a transit time of saidultrasound signal through said bone tissue.
 3. The method of claim 1wherein said processing step includes determining a velocity of saidultrasound signal through said bone tissue.
 4. The method of claim 1wherein said processing step includes determining an attenuation slopeassociated with said ultrasound signal.
 5. The method of claim 1 furthercomprising the step of independently directing said ultrasound signalfrom said first transducer to said second transducer through a mediumwith known acoustic properties and path length and free of said bonetissue to produce a reference electrical output signal and wherein saidprocessing step includes establishing a set of bone-oriented parametersassociated with said bone-oriented output signal and a set of referenceparameters associated with said reference signal and subjecting said setof bone-oriented parameters and said set of reference parameters tocomparative analysis.
 6. The method of claim 1 wherein said secondtransducer is said array transducer.
 7. A method of non-invasive andquantitative assessment of the status of bone tissue in vivo for atleast one of the quantities, bone-mineral density, bone strength, bonefracture risk, bone architecture and bone quality comprising the stepsof: positioning a first and second transducer distant from the skin onopposite sides of said bone tissue to provide gaps therebetween allowingfree insertion of said bone tissue unobstructed by the first and secondtransducers; acoustically coupling a pair of transducers to a compliantmedium bridging the gaps between the first and second transducers andthe distant skin; acoustically coupling the compliant medium to thedistant skin on opposite sides of said bone tissue; generating anultrasound signal and directing said ultrasound signal from onetransducer of said pair of transducers to another transducer of saidpair of transducers through said bone tissue, to produce a bone-orientedelectrical output signal; independently directing said ultrasound signalfrom said one transducer to said another transducer through a mediumwith known acoustic properties and path length and free of said bonetissue to produce a reference electrical output signal; establishing aset of bone-oriented parameters associated with said bone-orientedoutput signal and a set of reference parameters associated with saidreference signal; and subjecting said set of bone-oriented parametersand said set of reference parameters to comparative analysis, whereby anestimate of said at least one of the quantities, bone-mineral density,bone strength, bone fracture risk, bone architecture and bone quality isobtained.
 8. The method of claim 7 wherein said another transducer is anarray transducer.
 9. The method of claim 8 wherein said one transduceris an array transducer.
 10. A method of non-invasive and quantitativeassessment of the status of bone tissue in vivo for at least one of thequantities, bone-mineral density, bone strength, bone fracture risk,bone architecture and bone quality comprising the steps of: (a)positioning a first and second transducer distant from the skin onopposite sides of said bone tissue to provide gaps therebetween allowingfree insertion of said bone tissue unobstructed by the first and secondtransducers; (b) acoustically coupling a pair of transducers to acompliant medium bridging the gaps between the first and secondtransducers and the distant skin; (c) acoustically coupling thecompliant medium to the distant skin on opposite sides of said bonetissue; (d) generating an ultrasound signal and directing saidultrasound signal from one transducer of said pair of transducers toanother transducer of said pair of transducers through said bone tissue,to produce a bone-oriented electrical output signal; (e) repeating saidstep (d) a plurality of times to obtain a plurality of bone-orientedoutput signals; (f) averaging said plurality of bone-oriented outputsignals to obtain an averaged bone-oriented output signal; (g)independently directing said ultrasound signal from said one transducerto said another transducer through a medium with known acousticproperties and path length and free of said bone tissue to produce areference electrical output signal; (h) repeating said step (g) aplurality of times to obtain a plurality of reference signals; (i)averaging said plurality of reference signals to obtain an averagedreference signal; (j) establishing a set of bone-oriented parametersassociated with said averaged bone-oriented output signal and a set ofreference parameters associated with said averaged reference signal; and(k) subjecting said set of bone-oriented parameters and said set ofreference parameters to comparative analysis whereby an estimate of saidat least one of the quantities, bone-mineral density, bone strength,bone fracture risk, bone architecture and bone quality is obtained. 11.The method of claim 10 wherein said another transducer is an arraytransducer.
 12. The method of claim 11 wherein said one transducer is anarray transducer.
 13. A method of non-invasive and quantitativeassessment of the status of bone tissue in vivo for at least one of thequantities, bone-mineral density, bone strength, bone fracture risk,bone architecture and bone quality comprising the steps of: (a)positioning a first and second transducer distant from the skin onopposite sides of said bone tissue to provide gaps therebetween allowingfree insertion of said bone tissue unobstructed by the first and secondtransducers; (b) acoustically coupling a first transducer and a secondtransducer to nearby skin on opposite sides of said bone tissue, whereinsaid second transducer is an array transducer; (c) generating anultrasound signal and directing said ultrasound signal from said firsttransducer to an element of said second transducer through said bonetissue to produce a bone-oriented electrical output signal; (d)independently directing said ultrasound signal from said firsttransducer to said element of said second transducer through a mediumwith known acoustic properties and path length and free of said bonetissue to produce a reference electrical output signal; (e) establishinga set of bone-oriented parameters associated with said bone-orientedoutput signal and a set of reference parameters associated with saidreference signal; (f) repeating said step (c), said step (d), and saidstep (e) at least one time wherein said ultrasound signal is directed insaid step (c) and said step (d) from said first transducer to apreviously unselected element of said second transducer to therebycreate la plurality of sets Q_(s) of bone-oriented parameters and aplurality of sets of reference parameters; (g) evaluating said pluralityof sets of bone-oriented parameters to locate an anatomical region; and(h) subjecting at least one set of said plurality of sets ofbone-oriented parameters, said one set corresponding to said anatomicalregion, and at least one set of said plurality of sets of referenceparameters to comparative analysis, whereby an estimate of said at leastone of the quantities, bone-mineral density, bone strength, bonefracture risk, bone architecture and bone quality is obtained.